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Abstract 

We describe a novel artificial dielectric material which has applications at millimetre 
^ ■ and submillimetre wavelengths. The material is manufactured from layers of metal mesh 

Oj, patterned onto thin polypropylene sheets which are then bonded together using a hot 

JL I pressing process to provide planar rugged discs which can be reliably cycled to cryogenic 

in ' temperatures. The refractive index of this material can be tuned by adjusting the geom- 

c/2 . etry and spacing of the metal-mesh layers. We demonstrate its usage by designing and 

characterising a broadband anti-reflection coating for a Z-cut crystalline Quartz plate. 
The coating was fabricated and applied to the quartz using the hot press technique and 
V^ . characterized using a Fourier Transform Spectrometer. The performance is shown to be 

CD ' in good agreement with HFSS and transmission line modelling results. 
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1 Introduction 



Metal mesh technology has been adopted in the millimetre and sub-millimetre wavelength range 
(THz frequencies) as the standard method for fabricating optical filters, beam-splitters and 
dichroics [H |21 El H] . Components fabricated using this technique have been used in a number of 
/\ THz cameras and spectrometers El El El E] • More recently, the same technique has been used 
(^ to manufacture half-wave plates using aligned metallic patterns with air-gap spacing between 



the layers [TO]. These traditional metal mesh components are not considered metamaterials 
because their electromagnetic properties are not independent of their thickness, i.e. they cannot 
be characterized as having a bulk electric permittivity and magnetic permeability. 

This paper demonstrates that closely spaced multiple layers of metal-mesh films embedded 
in polypropylene behave as an artificial dielectric metamaterial. By changing the metal mesh 
geometry the effective impedance of the mesh layer is modified changing the apparent refractive 
index of the overall material. In Section 2 we present a theoretical analysis to obtain a physical 
insight into this behaviour and to show how the geometrical parameters modify the effective 
dielectric constant. This modelling uses a commercial finite-element analysis package: High 
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Frequency Structure Simulator (HFSS; [H]), which provides an exact solution to Maxwell's 
equations for a unit cell of the structure. 

As a demonstration, a prototype broadband anti-reflection coating was designed and applied 
to a single Z-cut quartz plate. In Section 3 spectral measurements of both the anti-reflection 
coating and of the quartz plate before and after coating with the artificial dielectric are pre- 
sented. These results are compared with predictions from HFSS and our own transmission line 
model. 

2 Theory and modeling 

2.1 Theoretical elements 

The metal mesh technology described in Ade et al. |lj is based on standard two port microwave 
circuit analysis. The metal mesh layer is modelled as an equivalent lumped LC filter and the 
spacing between meshes is modelled as a transmission line with impedance, Z = Zo/n where 
Zq is the impedance of free space and n is the index of refraction of the material between 
meshes. The lumped element impedance for a single sheet containing periodic metal structures 
has been previously determined by a number of authors [121 [131 [IH [IS] to be dependent only 
on the geometric properties of the metal. 

For filters and other components built up of many individual metal meshes, the metallized 
layers are stacked together with plane parallel spacers, which consist of either air gaps or 
dielectric discs. The model is successful in designing and predicting the performance when the 
spacing between layers d, is large enough so that there is relatively little inductive or capacitive 
coupling between layers. Air gap devices require an annular support ring, while the dielectric 
spacers and metalhzed sheets can be fused together using a hot pressing technique to make a 
solid self-supporting disc. While the design of these components is based on radiation at normal 
incidence to the metal mesh layers, their performance, which is measured in a converging beam 
of half angle 9.6 degrees (see later), agree well. This indicates that there is not a strong 
dependence with off- axis angle. 

The simplest metal mesh pattern used for low pass filters is termed a capacitive grid as 
described in Ade et al. [1]. Here we describe the properties of a stack of identical capacitive 
mesh grids with spacing between the grids, d << X immersed in a dielectric substrate layer 
using hot-pressed technology. In these structures, the two port microwave circuit model breaks 
down due to the capacitive coupling between layers. However, the properties of the components 
are well described by an effective bulk permittivity, €,.(//) > 1 and bulk permeability, /ij. = 1 
corresponding to a material with effective index of refraction, neg > 1. By changing the 
basic mesh geometry, the mesh layer spacing and the number of mesh layers the effective 
permittivity of the newly created artificial dielectric slab can be modified. This is expected 
since the lumped capacitance of the grids will increase the capacitance per unit length for an 
electric field parallel to grids, thus increasing the effective permittivity of the material. A full 
EM analysis of structures made of solid metal bars with gaps between them has been reported 
[T6] . We perform a semi-empirical analysis of our structures based on 3-D EM simulations using 
HFSS and compare these results to measurements of fabricated components. 



2.2 Essential parameters in the build 

Figure 1 shows the basic structure of the capacitive metal mesh layers embedded in a base 
dielectric material such as polypropylene. 
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Figure 1: The basic structure of the artificial dielectric. The structure shown has all of the 
metal meshes in the stack aligned with each other. Small offsets in alignment from one layer 
to another have a small effect on the effective index of refraction. 

The structure of a single mesh layer is periodic so there are only three essential geometrical 
parameters in the model; the spacing between the layers, d, the gap between two adjacent 
metallic squares, 2a, and the period of the repetitive square structure, g, as shown in Figure 2. 
It was also found useful to categorize the mesh patterns by the ratio a/g since for a given grid 
period this characterizes the basic frequency dependence of the single mesh. 



2.3 HFSS simulation result for metal mesh dielectric 

We used HFSS simulations to explore the optical properties of these capacitive grid stacks with 
different geometrical properties. For these simulations we fixed the basic grid periodicity, g, 
to be 100 /im. Changing this dimension along with the spacing between the layers scales the 
observed properties to different spectral regions. The simulations were made for a range of a/g 
ratios from 0.05 up to 0.28 (0.05, 0.1, 0.14, 0.28). For each a/g, we also varied the number of 
layers (5, 10, 20 layers) and explored the dependence on the grid spacing, d, between adjacent 
layers (4, 5, 8, 10, 12, 15 and 20 /xm). 

The simulation transmission results for one set of a ten layer device with a/g = 0.28 for 
various grid spacing are shown in Figure 3. 

The refractive index, n, is derived from the transmission data by assuming that the resultant 
material behaves like a plane parallel dielectric plate with thickness. Ax. From the HFSS 
simulation we obtain the complex reflection and transmission coefficients, r and t, of the material 
and invert the standard Fabry-Perot (FP) transmission formulas [T7] : 
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Figure 2: Top view of single layer metal mesh pattern. The evaporated Cu squares created 
photolithographically are typically 0.1 /im thick. 
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Figure 3: HFSS simulated transmission for 10 mesh layers immersed in polypropylene with fixed 
grid geometry a/g = 0.28 and g = 100 ;um for mesh layer spacings of 4 (rightmost minimum), 
8, 10, 12, 15 and 20 /im (leftmost minimum). 
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to determine the index of refraction, n and material impedance, z. From these we can 
compute values for the relative permittivity, e(o") and permeability, ii[(y\ We find for all 
structures that /i(cr) ~ 1 independent of frequency so the material can be characterized as a 



simple dielectric and the index can be calculated from the FP intensity formula 

1 
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1 + [2r/(l - r2)]2sin2(27rn(TAa;) 
or simply by measuring the transmitted intensity at the first minimum: 
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Inversion of this for the transmission at the first minima observed in Figure 3 gives the 
effective refractive index for the artificial material at that frequency. 
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Figure 4: Predicted refractive index n for a 10 mesh layer device as a function of mesh layer 
spacing, for mesh geometries a/g of 0.05, 0.1 and 0.28, at 5cm"^ wavenumbers (frequency 
150 GHz), with g = 100 /xm. Errors are 2% due to simulation accuracy. 

The estimated refractive index as a function of the grid spacing for each a/g ratio is presented 
in Figure 4. This shows that for a fixed a/g ratio, the refractive index of the material decreases 
as the grid spacing increases. It can also be seen that for a fixed grid spacing, a higher a/g 
ratio yields a lower refractive index. For both cases an increase in the grid spacing or the a/g 
ratio results in a decrease per unit length of the effective artificial material capacitance for an 
E-field parallel to grids and hence a decrease in refractive index. We perform simulations of 
structures with different numbers of mesh layers and find that the effective index of refraction is 
independent of component thickness and therefore can be considered a bulk material parameter. 
Finally, Figure 5 shows that the index of refraction is weakly dependent on the frequency and 
that the metamaterial behaves as an artificial dielectric over a wide range of wavelengths, 
corresponding to 5^ < A. 



2.4 HFSS model parameters for anti-reflection coating 

One important application of this artificial dielectric is as an element in a broadband multilayer 
anti- reflection coating (ARCs). In this application there are two fundamental material require- 
ments that need to be met. First, the materials used in ARCs must have a range of appropriate 
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Figure 5: Predicted refractive index n for 10 mesh layer device and fixed spacing 10 /im, as a 
function of frequency, for mesh geometries a/g = 0.05, a/g = 0.1, a/g = 0.28, with g = 100 /xm. 
In general, there is a slight increase in the effective index of refraction with frequency due to 
the increase in spacing between the layers relative to the wavelength. 



refractive indices and high transparency over the desired spectral band. Secondly, the material 
must be mechanically suitable for applying to crystalline materials. 

ARCs are usually used to maximize the device transmission over spectral bandwidths ap- 
proaching 100 %. Previous designs used specially prepared polypropylene layers p^ [20] that 
are loaded with high refractive index powders (Ti02) or ceramic based materials (e.g., Rogers 
TMM material [21]) to create a particular refractive index. These ARCs needed three layers 
to achieve a flat response across the band. Each layer required hot pressing or grinding to the 
required thickness before being bonded to the front and back surfaces of the crystalline plates 
fi9[ [22] . Major problems were encountered with both the previous approaches. The loaded 
powder layers were slow to manufacture since the powder needed to be uniformly mixed in the 
polypropylene and then the layers had to be hot pressed to the appropriate thickness. The 
ceramics are brittle and could only be thinned by using a grinding technique which proved to 
be time consuming and difficult to handle in thin sheet form. 

Compared with the previous techniques, the artificial dielectric design has the advantages 
of simple refractive index tuning through the geometry and spacing of the grids. Without 
high index powders added into the material, the metal mesh design has a complete control 
over the thickness of the coating layer and the coating material is not brittle and has better 
performance in thermal cycling. Furthermore, since the base dielectric is polypropylene, it also 
has low absorption for frequencies up to 20 THz [23]. Thus the artificial dielectric meets with 
the ARC fundamental requirements. 

To demonstrate the new material we have prototyped an ARC for a Z-cut crystal quartz 
plate with refractive index n=2.1. To increase the spectral bandwidth of the low refiection 
region we needed two materials with intermediate refractive index values close to 1.3 and 1.7. 
The first of these is realised using readily available porous PTFE whilst the latter requires the 
customised artificial dielectric. 



The artificial dielectric is constituted by two essential elements, a metal mesh grid and 
the embedding dielectric. The first is obtained by photo-lithographic techniques adopted in 
the manufacturing of far-infrared low-pass filters |24j . The copper deposited patterns on thin 
substrates (8/xm) of polypropylene are then inserted into a stack (see Figure 6). The two metal 
mesh layers were immersed in the substrate at a distance of 8/xm from the top and bottom 
surfaces and with a spacing s=24/im between two layers. The single mesh layers are periodic 
structures of square grids patterns. In this design a/g = 0.14 with g=25Afim. and a=3.556/im. 
Once the 40 /im multi-layer structure is assembled it is then fused by hot-pressing the layers 
at high temperature near the polypropylene melting point (160° C). The temperature is such 
that the polypropylene layers are joined but is insufficient to allow the polymer chains to flow 
in a liquid-like manner. Expansion due to relaxation when hot pressing the layers is avoided 
by maintaining high pressure in the press and through the absence of soft materials in the 
stack. Thicknesses measurements made after hot pressing show that the layer thicknesses are 
maintained to better than ±1 /xm. 



Metal Grid 
Substrate 




Figure 6: Prototype ARC metal mesh model. 

The artificial material created by this process has the consistency of a solid plastic film 
which can be easily handled and cut to size to suit requirements. Bonding of this material to 
the quartz crystal is again performed through a cycle of thermal heating with applied pressure 
using a thin layer of polymer based glue (thickness << wavelength). Finally, the outer layers of 
porous PTFE are again glued at the extremes of the stack (top layer in Figure 7)with identical 
procedure (the two procedures can in fact be performed in one single press cycle). 

The ARC layers are symmetrical with respect to the central crystal in order to taper the 
impedance that a given wavelength of radiation experiences when entering and subsequently 
exiting the crystal. 



2.5 HFSS prediction for anti-reflection coating 

Simulations for both the prototype ARC on its own and the completed quartz plate with AR 
coatings have been made. The transmission for the prototype ARC layer is shown in Figure 8. 
The refractive index was calculated based on transmission as described in section 2.3 to be 1.65 
in average. Figure 9 shows the comparison between the uncoated quartz substrate, the AR- 
coated substrate and the simulated transmission expected. The efficiency of the coating yielding 
high transmission over the frequency band from 18.6-52 cm~^ (560-1560 GHz) is obvious. 
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Figure 7: Full model ARC on both sides of quartz substrate. 
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Figure 8: Transmission for ARC metal mesh layer: measurements (black dotted line), HFSS 
simulation (black line with diamond symbols) and Fabry-Perot model (grey line). 



3 Measurements and Data Analysis 

As a stand-alone material, the embedded grids in polypropylene will have a complex impedance 
with the real part mimicking a behaviour analogous to that of a plane parallel dielectric slab. 
We therefore subjected the material to the same spectral tests that we would have imposed on 
a thin dielectric slab to determine its refractive and absorptive properties. 

The material has also been measured when used as a tuned anti-refiection coating element. 
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Figure 9: Transmission of the complete ARC quartz plate. Lower grey line is measured data 
of the uncoated quartz substrate, the black line is the HFSS simulation of the complete ARC 
quartz plate and the black dotted line is the measured data of the complete ARC quartz plate. 

As reported above for this measurement it was hot pressed on both sides of a quartz plate. We 
were careful to choose grids with a period of 25.4 yum to ensure that diffraction effects would be 
negligible at the frequencies for which the ARC was optimised. However, we expected to ob- 
serve spectral features resulting from the complex impedance nature of the resonant structures 
which would distinguish the measured spectral behaviour from that of a normal dielectric with 
constant index of refraction. 



3.1 Spectral measurements of the material. 

To verify the optical properties of the artificial dielectric alone and with the crystal, we placed 
the material in the //3 beam of a rapid scanning polarizing Fourier Transform Spectrometer 
(FTS). The spectral features associated with the source (a Hg-arc lamp) together with the 
spectral efficiency of the detector (a 4 Kelvin Ge bolometer) were removed by performing the 
ratio with a background spectrum recorded with the sample removed. Three sets of 10 scans 
each were performed in order to obtain acceptable signal to noise for a spectral resolution 
of 0.25 cm^^ and limit phase errors in the Fourier transform process. Interferograms in each 
set were averaged and then Fourier Transformed. The resulting three spectra have then been 
averaged before ratioing against the background spectrum. The final transmission spectrum is 
shown in Figure 8. By repeating this measurement at various orientation of the mesh axises with 
respect to the polarization direction we confirmed there is no polarization sensitivity observable 
above 1 % level as expected. 

From the peak of transmission at 62cm~^ we can infer that absorption loss due to the 
polypropylene, which is usually a monotonic function of frequency, is negligible. Data for the 
absorption coefficient of Tucker and Ade [23] indicate that the expected loss would be only 
0.2 %. With this consideration we can expect the transmission of an equivalent dielectric slab 
to be given by equation (1). By inverting this equation whilst fixing the physical thickness t 
of the material we can numerically solve for each value of a to obtain a value n{a). Figure 



10 shows the residual mismatch of the FTS data with the ideal Fabry-Perot behaviour. The 
contour levels are the percentage error on the transmission mismatch (measured/ideal). The 
thick line in the clear contour region shows the best fit over the spectral interval 20-40 cm^^ 
for an assumed linear function of frequency. The linear fit parameters are: 



n = uq + uiu 
no = 1.637 
ni = 2.07 X 10"''^ 



(5) 
(6) 
(7) 
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Figure 10: Residual contour of measured transmission vs expected ideal behaviour of an ideal 
dielectric slab. The thick line in the clear region denotes the best fit linear function of frequency. 

The small ripple artifacts which are visible in the leftmost region of the contour, arise from 
standing waves in the FTS sample box between the sample and the output polariser, and hence 
have been ignored in the fitting procedure. 

Using these values as initial conditions for a simulation of the complete stack we can then 
proceed to verify that its employment as an anti-refiection coating element is consistent with 
its stand-alone properties. 

3.2 Spectral measurements of coated quartz substrate. 

As previously stated, the end goal is to produce a broadband anti-reflection coating to remove 
reflection loss from high refractive index crystalline components. The complete ARC structure 
as described above was measured in the identical setup as the stand-alone material. Figure 9 
presents these data which can now be compared with the spectral performance as simulated at 
a number of frequencies using HFSS. 

First, we note that in general there is good agreement between the envelopes of the fringes. 
This denotes that the material is acting as a good anti-reflection coating over the desired band. 
However, some minor discrepancies are also present which question either the set of input 
parameters for the simulation or our understanding of the effects. 
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Having manufactured a large number of anti-reflection coated windows and filters with 
different techniques, we are aware that some of the parameters can slightly change due to 
the heat-bonding process in our press. Indeed some bonding materials change when heated. 
Porous PTFE can be pressed to a smaller thickness and being porous can absorb some of the 
polyethylene glue, raising the equivalent index of refraction by a small amount and changing 
its final thickness. Polypropylene on the other hand tends to relax and expand if there is not 
sufficient pressure placed on it. Crystals clearly will not be deformed at these temperatures 
and pressures. 

Performing HFSS simulations with different values of parameters at multiple frequencies 
proves to be very time consuming, so we have adopted a transmission line based model [22] to 
perform a Monte-Carlo simulation cycling through some of the input parameters (thickness of 
the coating layers and indexes of refraction to find the solution which closely matches our data 
set and check for potential manufacturing effects. 

In order to justify the results we first tested the consistency of the two simulation/modelling 
packages by running the transmission line model on the nominal input-parameters. In Figure 
11 we can clearly see an excellent agreement from low to high-frequencies with the two best- 
matching regions in 20-25 and 42-46 cm~^. The loss of relative precision in the HFSS simulation 
at the higher frequencies is potentially the cause of the mismatch occurring there. 




30 40 
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Figure 11: Expected match of the scattering matrix model (grey) with the HFSS simulation 
(black dotted). 

The Monte- Carlo parameter fit results and minimum square residuals are given in Table 1 
for the best fit of the data shown in Figure 12. This fit, performed over the region 15-55 cm~^, 
yields a good match with the following caveats regarding the numbers in Table 1: 

• Mupor [25]: Porous PTFE. The increase in the mupor index and its decrease in average 
thickness is consistent with compression of a porous material during bonding. 

• ADM - t: The thickness of the Artificial Dielectric Material (ADM). The values of thick- 
ness slightly larger than expected could point to thermal relaxation of the polypropylene 
layers. This phenomenon has been observed in direct coating of quartz windows. This 
tendency of thin layers to shrink in diameter and grow in thickness of a few percent can 
occur when the pressure applied is not sufficient or other less dense materials in the stack 
give way. 
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Figure 12: Transmission spectrum of the stack (PPTFE-ADM-Quartz-ADM-PPTFE). The 
black hne is the best fit of a scattering matrix model with varying optical constants to fit the 
behavior of the ADM as a function of frequency. 

• ADM - n: The refractive index of the Artificial Dielectric Material (ADM). The values 
of index of refraction are slightly different from the recovered functional values found 
previously (no = 1.62 and rii = 2 x 10~^) but it is important to note that in the interval 
of frequencies considered, the two functional dependencies are not too dissimilar. 

• Top-bottom asymmetry : Generally we would not expect a marked asymmetry between 
top and bottom coating layer thicknesses. However, we note that in the vacuum oven the 
top plate of the press is only conductively heated through the stainless steel press screws. 
A temperature difference of a few degrees will affect the tendency of materials to change 
their thickness subject to pressure. 

• Quartz n : The values of quartz index obtained are slightly different (smaller) than those in 
literature [26] by a scaling factor ~ 1/1.007. This is expected as the standing wave pattern 
in the fitted data does not vary with a constant product of n and the quartz thickness. To 
avoid additional standing waves in the measuring setup, the sample was tilted by i =~ 
10° ±5° and this would affect the thickness travelled by l/cos[sin~^(n • sinz)] ~ 1/(1 + e) 
with e C [1,8] ■ 10^^. By forcing the nominal thickness the fitted value of the index of 
refraction compensates accordingly. 



4 Conclusions 

We have shown through modelling and experimental verification that artificial dielectrics with 
refractive indexes above that of the base material can be obtained over a broad spectral band 
by fusing in layers of metal mesh. The resultant refractive index can be easily controlled by 
adjusting the geometrical parameters of the meshes and the spacing between them in the new 
material. 

We have also demonstrated the viability of this metamaterial as an element in producing 
a broad band anti-refiection coating for a quartz crystal. Again modelling and measurement 
results are in good agreement showing that we have good control of this process. As a further 
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Table 1: Simulation results. 

List of parameters (fixed and variable) in the simulation of the coated stack. It is important 
to specify that due to the nature of the fit and the different parameter intervals, the errors 
reported are only indicative of the dependence of the result from that parameter (i.e. how 
much does that variable influence the result) and not a simple standard deviation of a gaussian 
distribution. The two errors are respectively indicating an increase of (1)=10 % and (2)= 100% 
in the residual of the least square quadratic flt when applied to the speciflc parameter. (*) The 
quartz value was measured on the substrate with a precision micrometer and the error resulting 
in the repeatability of the measurement across the sample. 



Parameter 


# 


Value 


Error^i^i) 


Error(2)(±) 


Units 


Mu-por Index of Refraction 


^fipar 


L375 


0.01 


0.03 




Mu-por Thickness (Top) 


fipor 


4L9 


1.5 


4 


(/im) 


Mu-por Thickness (Bottom) 


lipor 


44.1 


1.7 


4 


(^m) 


ADM Index of Refraction 


no 


1.592 


0.014 


0.04 




ADM Index of Refraction 


rii 


7.5e - 3 


5e-4 


1.4e - 3 




ADM Thickness (Top) 


'^ADM 


42 


+2 
-1 


+5 
-2 


(^m) 


ADM Thickness (Bottom) 


^ADM 


45 


+2 
-1 


5 


(/im) 


Quartz Index of Refraction 


no 


2.095 


0.002 


0.005 




Quartz Index of Refraction 


rii 


4e-6 


4e-5 


1.4e-4 




Quartz Index of Refraction 


n2 


2e-7 


le-6 


4e-6 




Quartz Loss Tangent 


«! 


4.1e-5 


7e-6 


2e-5 




Quartz Loss Tangent 


^2 


4.0e-6 


2e-7 


6e-7 




Quartz Thickness (*) 


(i(flxed) 


2.890 


— 


0.003 


{mm) 



test of the robustness of the coated layers we dropped the sample into liquid nitrogen (77K) 
and warmed to room temperatures flve times with no discernable signs of the coating peeling 
of fracturing. 

We believe that this new approach will assist in the manufacture of many other optical 
components in the FIR to millimetre wavelength range. These include flat lenses, birefringent 
components for waveplates (using 1-D grid structures) and dispersive components. 
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